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ON CELL ENTROPY INEQUALITY FOR DISCONTINUOUS 

GALERKIN METHODS 1 


Guangshan Jiang and Chi- Wang Shu 
Division of Applied Mathematics 
Brown University 
Providence, Rhode Island 02912 


ABSTRACT 

We prove a cell entropy inequality for a class of high order discontinuous Cialerkin finite 
element methods approximating conservation laws, which implies convergence for the one 
dimensional scalar convex case. 
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Grant NAG-1-1145 and Contract NAS1-19480 while the second author was in residence at ICASE, NASA 
Langley Research Center, Hampton, VA 23681-0001, and AFOSR Grant 93-0090. 
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1 Introduction 


In [3] Cockburn and Shu defined a class of discontinuous Galerkin finite element methods 
for conservation laws (the multi-dimensional case was discussed in [4]): 

u t + div f(u) = 0 (1-1) 

The entropy solution of (1.1) also satisfies 

U{u) t + div F{u) < 0 (1-2) 

in the distribution sense, for any convex function U{u) and consistent entropy flux F(u) sat- 
isfying F'(u ) = U'{u)f'{xi). The scheme is obtained in the following way: first a triangulation 
T h is chosen. In one space dimension is just a collection of subintervals Ij = (:Cj_i, £ J+ |), 
which are not necessarily of the same length. We will use the notation A x 3 = Xj+l ~ 
h = maxj Ax.j and x : = 1(^1 + x j+ i). In two and more space dimensions T h is a collection 
of triangles, tetrahedrons or other simple geometric objects. The solution space 14 is defined 
as the collection of all piecewise polynomials of degree up to r for a (r + l)-th order method. 
The functions in V/, are allowed to have discontinuities across element interface. The conser- 
vation law (1.1) is multiplied by a test function v £ 14, integrated over an element I< € T /t , 
and formally integrated by parts to shift the spatial derivatives from f(u) to v. The result 

is: 


/ u t (x,t)v(x)dx+ I f{u{x,t))-nv(x)ds - f{u(x, t)) ■ grad n(x)dx = 0. (1.3) 

JK e€dJ< Je JK 

where dl< is the boundary of the element I< and n is the unit outward normal vector. This 

procedure up to now is just the standard discontinuous Galerkin method. Two features from 

the finite difference methodology are then used in [3], [4]: one is the monotone or E flux 

(approximate Riemann solver for systems) to define the trace of the physical flux f(u) at 

the cell interface; the other is a local nonlinear limiter which limits the values of u at the 
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cel] interface by the differences of the means of u over elements. With these two ingredients, 
it was proven in [3] and [4] that the schemes are formally high order accurate, measured by 
local truncation errors in smooth regions including at extrema, are total variation bounded 
for one space dimension, and are maximum norm bounded for any space dimensions. The 
only thing missing is the entropy condition: in [3], we were able to prove entropy consistency 
for the square entropy U(u) = ^ for one dimensional convex f(u) with a ^-independent 
modification to the scheme, following the idea of Osher [12]. We were also able to prove 
entropy consistency for all convex entropies with a ^-dependent modification to the scheme. 
The /i-dependent limiters make the proof of high order schemes easy, but it is not very desir- 
able for practical computations, because it usually limits the slope near the discontinuities 
stronger than necessary, and essentially flattens the solution there to piecewise constants for 
fine mesh. It also destroys the self-similarity of the scheme. 

Entropy condition seems difficult to prove for high order finite difference schemes. Osher 
and Tadmor proved [13] that finite difference schemes (which evolve only the mean) which 
satisfy cell entropy inequalities for all convex entropies can be at most first order accurate. 
Even for one entropy inequality (say for the square entropy), for one space dimension and 
for convex /(«), the proof is extremely elusive if one does not modify the scheme. Osher [11] 
and Nessyahu and Tadmor [10] were able to get such cell entropy inequalities, for the square 
entropy and for the second order MUSCL scheme, with some ^-independent modifications; 
Yang [14] was able to prove convergence of unmodulated second order MUSCL scheme 
using a global analysis rather than relying solely on cell entropy inequalities; Lions and 
Souganidis [8] proved convergence of second order MUSCL scheme for steady state Hamilton- 
Jacobi equations and conservation laws. There are also many results which prove entropy 
consistency and/or convergence using /i-dependent limiters or modifications for high order 
schemes. For example Coquel and LeFloch [7] for finite difference; Johnson, Szepessy and 
Hansbo [9] for streamline diffusion finite element; Cockburn, Coquel and LeFloch [5] and 
Cockburn and Gremaud [6] for high order finite volume, streamline diffusion or discontinuous 
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Galerkin schemes with /i-dependent “shock capturing” terms. These results are usually more 
general (multi space dimensions, non-convex fluxes, etc.). However, as we have indicated 
before, /(-dependent limiters or modifications should be avoided for practical computations if 
possible. More recently, Bouchut and Perthame [2] obtained a second order one dimensional 
scheme which is consistent with all entropy conditions and does not use /i-dependent limiter. 
A key ingredient of [2] is to evolve both the mean and the slope and use the whole function 
(not just the mean) to obtain cell entropy inequalities. Discontinuous Galerkin methods also 
fall into this category (evolving the whole polynomial in a cell rather than just the mean), 
and we use this fact strongly in this paper. 
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In Section 2 we prove a cell entropy inequality using the square entropy U(u) — y 
for the unmodulated semi-discrete discontinuous Galerkin method of [3], [4]. The proof is 
remarkably simple and does not even use any nonlinear limiters. It works thus for any spatial 
order of accuracy. As far as we know, all the previous cell entropy inequalities without h- 
dependent limiters must restrict the slope of a function by £^mmmod(A+Uj, A_u ; ) (which 
means that the difference at the interface, ut , j — u i_ , must be of the same sign as that of 

J T 2 j' 2 

Uj +1 - Uj, i.e., no “sawtooth” is allowed in the reconstruction), hence cannot be higher than 
second order accurate. The result in this paper illustrates the potential of discontinuous 
Galerkin methods, or equivalent Hermite type finite difference/finite volume type methods, 
which evolve the whole polynomial in the cell rather than just the mean. Time discretization 
is discussed in Section 3. 

2 Cell Entropy Inequality for the Square Entropy 

The discontinuous Galerkin scheme in one space dimension, defined in [3], in its semi-discrete 
form without slope limiting, is the following: Find u(-,f) 6 Vh, such that, for all v E U and 
all subinterval /j: 


I itj(x, t)v(x)dx + hj + i(t)v(Xj + i) — hj_i(t)v(x+_i) — f(u(x,t))v x (x)dx — 0 (2.1) 

J i j ^ ( 
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Here, h- + i(t) = h (u{x~ + i,t),u(x+ +l ,t)J is a Lipschitz continuous monotone flux (i 
non-decreasing in the first argument and non-increasing in the second argument), o 
generally an .E-flux as defined by Osher [11]: 

u + ) — /(u)) ( u + — u ~ ) < 0 

for all u between u~ and u + . Some examples of the commonly used monotone fluxes 
found in, e.g., [3]. 

If we take u(ar) = u(x,t ) in (2.1), we get 


/.PH 1 ) 

We define 


dx + h J+ i(i)u(x j+L ,t) - h j _i(t)u(x+_ 1 ,t) - j f(u(x,t))u x (x,t)dx = ( 

2 2 J J j 


9(u) = J f(u)du 


and rewrite (2.3) as 


where 



dx + F j+ i(t) - E,_i(<) + Aj{t) = 0 




t) - g{u(x j+ ^t)) 


is consistent with the entropy flux for the square entropy: 


and 


F(u) = J f'(u)udu = f{u)u - J f(u)du - f(u)u - </(u), 


Aj(t) - —hj_i(t) (u(x+_±,t) - u(x^_i,f)^ +if(u(xt_i,<)) - g(u(x._±,t)) 


.e. h is 
r more 

( 2 . 2 ) 
can be 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 
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(2.8) 


= -(h H m-m) o) >o 

where we have used the mean value theorem and the definition (2.4) of g(u) in the second 
equality, £ is between and t), and the last inequality is due to the property 

of E-flux (2.2). 

We have thus proven the cell entropy inequality 


l dx + P i+ L(t) - F H (t) < 0 (2.9) 

for the square entropy U(u) = f . Notice that we do not need any nonlinear limiting at 
this stage. However, nonlinear limiting as introduced in [3] and [4] will not destroy this cell 
entropy inequality (see next section). The cell entropy inequality (2.9) trivially implies L 
stability of the scheme (again without even using the nonlinear limiting). 


d_ 

dt 



dx < 0 


( 2 . 10 ) 


but it is much stronger. For example, if f(u) is convex and we use the nonlinear limiting [3] 
to obtain a total variation boundedness for the solution, we will have convergence towards 
the unique entropy solution. 

The same entropy inequality can be obtained for multi-space dimensions with arbitrary 
triangulations: 


I dx+Y. [F.M*,t)ds< 0 

Jk \ 2 ) t eedK Je 


( 2 . 11 ) 


where F e ,/c is consistent with F ■ n K for the entropy flux F in (1.2) and the outward normal 
n K , and F e k — —F ey K‘ f° r th e f wo neighboring elements K fl^ — l e ) (conservation). We 
omit the derivation detail since it parallels that for the one dimensional case. 


We have the following two remarks: 

(1). If we try to do the same estimate for a general convex entropy U(u), we can obtain 
exactly the same cell entropy inequality module an interpolation error term: 


5 



- / («(«, 0, + /(«(*, 0)*)(^(«(a:,<)) - n U'(u(x,t)))dx (2.12) 

where Ilu is a projection into the space V), which interpolates at the two end points of Ij. 
This motivates the following ^-dependent modification to the scheme, which is similar to the 
“shock capturing term” added to the streamline diffusion method in Johnson, Szepessy and 
Hansbo [9] and in Coc.kburn and Gremaud [6]: 


J^u,(x,t)v(x)dx + h j+i (t)v(x- +j ) - VlW»(*, + _i) - l f(u(x,t))v,(x)dx 

+ Chj M*, 0 + /(«(*, 0W rr^ = 0 ( 2 . 13 ) 

‘3 yii x yx ? t j I 

where C is a suitable positive constant. This, together with a L ^ bound which can be 
obtained by using nonlinear limiters [3], will give us a cell entropy inequality for arbitrary 
convex entropy at least for the r = 1 (second order) case, hence convergence for any noncon- 
vex flux f(u) in this case. For general r, more /i-dependent modification is needed. Notice 


that by adding this modification, the formal order of accuracy of the scheme is not changed: 
+ f(u(x,t)) x | is just the local truncation error. A similar argument as in [9] shows 
that the modification does not destroy convergence towards weak solutions (conservation). 
(2). The so-called discontinuous Galerkin method can also be recasted as a finite differ- 


ence scheme (Hermite type: where one evolves both the mean and the slope, maybe more). 
For example, the second order case is just the following scheme for the mean u 0 {t) and the 
slope Sj (t) where u(x,t ) = uj(t) + sj(t)(x - x/) in cell Ij (see [3]): 


dt Sj W = ~Ax? [/i+iW + /j-i(0] + ^ 3 / f (uj(t) + Sj{t)(x- Xj ))dx (2.14) 

where 

~ h (“i(0 + ~ — ^ +1 s j+ 1 (t (2.15) 
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We have thus in effect proved cell entropy inequalities for such high order Hermite type 
finite difference schemes without using the help of any nonlinear limiting. Of course, to get 
convergence one must use the nonlinear limiting to obtain L ^ and/or total variation bounds, 

3 Time Discretization 

We discretize (2.1) in time by the following class of methods: 


where 



- J^f(u n+e (x)^v x (x)d 


x — 0 


(3.1) 


u n+ '(x) = (I -0)u n (x)-Mu" +1 (*) 

- + i>,“” +, (-; +i )) (3.2) 

For 6 = 0 this is the Euler forward discretization; for 6 = 1 this is Euler backward; and 
for 6 = | this is Crank-Nicolson. 

If we take v(x) = u 1l + e (x), we obtain just as before 



where 


with g(u) defined by (2.4). We can rewrite (3.3) as 


(3.4) 


1 , (U " I(I) L c - K ^ u " (x)) ^ * o («> 
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Thus a sufficient condition to get the cell entropy inequality 


f (“"^(s)) 2 — ( uM (?) ) ? . (fa + F n+ f - F n+ i < 0 

Ji, 2A t J+ 2 J-j " 


(3.6) 


is just 6 > i.e. implicit schemes from Crank-Nicolson to Backward Euler. 

Up to now we have not considered the nonlinear minmod type limiters in [3] and [4]. 
These limiters will render the scheme total variation stable (in one space dimension) or L ^ 
stable (in multi space dimensions). We have to make sure that these limiters do not destroy 
the cell entropy inequality (3.6). This turns out to be a simple issue: suppose u n+1 (x) is the 
solution obtained from the scheme (3.1) without limiter, then projected to get the solution 
at time level n + 1 by some, minmod type limiter it n+1 (x) = Pu n+1 (x) (this is the procedure 
adopted in [3] and [4]). The above derivation for the unlimited scheme will give (3.6) for 


f ( un+1 ( J ) _ ) j — + fn+e _ pn+e < Q (3.7) 

Ji } 2 A t >+k M “ 

hence a sufficient condition to get the cell entropy inequality (3.6) for the limited (projected) 
solution u n+1 is to require the projection P to satisfy 


/ ( Pw(x)) 2 dx < / (w(x)) 2 dx (3.8) 

Jlj Jlj 

for all polynomial w(x) of degree up to r (recall that r + 1 is the order of the scheme), i.e. 
P does not increase the L 2 norm of any r-th order polynomial in cell l y Notice that this 
is the idea used in [2]. An easy way to ensure condition (3.8) is to write tn(;r) and Pw(x) 
as expansions of (scaled) Legendre polynomials pi(x) as in [3]: m(x) = Y^i=o a iPi{ x ) an ^ 
Pw(x) = Y?i=c>biPi{ x )- A simple sufficient condition to ensure (3.8) is now just \b[\ < |«/| 
for all /, which is easily checked to be correct for the P 1 case (r = 1) in [3] and also for the 
general P T case with similar minmod limiters on «/ (see also [1] for such limiters). 

It is also possible to get cell entropy inequalities for certain explicit time discretizations. 
For example, with leap-frog time discretization: 
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